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Abstract 

We investigate the possibility that neutrinos form superfluid-type conden- 
sates in background cosmological densities. Such condensates could give rise 
to small neutrino masses and splittings, as well as an important contribution, 
perhaps, to the cosmological constant. We discuss various channels in the 
context of the standard model. Many of these do not support a condensate, 
but some mixed-flavor channels do. We also suggest a new interaction, acting 
only among neutrinos, that could induce a neutrino Majorana mass of order 
1 eV. 
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For various reasons, but especially the recent evidence for neutrino oscillations and hence 
neutrino mass [1], considerable effort has been devoted of late to neutrinos and their prop- 
erties. Furthermore, another more speculative consideration has piqued interest in small 
neutrino masses: a number of cosmological observations has led to renewed interest in the 
long-discarded notion of a cosmological constant [2]. The preferred value of this constant 
leads to a mass scale that is quite similar to the mass for neutrinos inferred from the at- 
mospheric and solar neutrino observations, in the range of tenths of an eV to a few eV [1]. 
It seems to us that a natural framework in which to relate these seemingly disparate phe- 
nomena would be the formation of some type of neutrino condensate, the energy of which 
might provide the cosmological constant, while the neutrino masses could emerge from an 
expansion about the symmetry-breaking vacuum. 

Appropriate tools to investigate the possibility of neutrino pair condensation of the su- 
perconductor type, have been recently employed to explore the possibility of qq condensates 
in hadronic media at sufficiently high densities. This has been discussed in the context 
of a four-fermi approximation to QCD, whether induced by instanton effects [3] or by one- 
gluon exchange [4]. Similar effects have also been studied in the large- N approximation in 
a 2-dimensional model related to the Gross-Neveu model [5]. 

For our purposes, we look for pairing phenomena in the electroweak theory. Specifically, 
we suggest that condensates involving neutrinos may occur in cosmological situations where 
the relevant chemical potentials are non-zero. Of course, the chemical potentials of the 
various neutrinos in the universe are presumably quite tiny [6], reflecting as they do the 
difference in density between neutrinos and anti-neutrinos of a given species. We shall 
return to the estimation of the magnitude of the effects we discuss below; first we offer 
a purely theoretical discussion of the following problem: to identify, if any, the attractive 
channels in the electroweak theory that might permit neutrinos to condense. It should be 
noted that the earliest speculation on this subject of which we are aware [7], took place 
before the emergence of the standard model and the discovery of neutral currents, and the 
authors simply assumed the existence of an attractive channel. 

The hallmark of this type of pairing is that, once a Fermi surface forms, any attractive 
channel, no matter how weak, can produce a condensate; this is in contrast to the more 
familiar situation in hadronic physics of a chiral condensate, which typically requires the 
coupling to exceed a certain threshold value. 

The reason for this behavior of the pairing condensate can be understood in a variety 
of ways. If one looks at the gap equation in a mean- field approximation one sees that, as 
the gap tends to zero, a singularity develops at the Fermi surface; hence a non-zero gap is 
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necessary to prevent the formation of this singularity. 

Alternatively, in a renormalization-group approach [8] one sees that the renormalized 
coupling tends to infinity as one integrates out all modes above and below the Fermi surface; 
thus even an apparently weak coupling grows without bound and allows the condensate to 
form. 

Our investigation will be limited to the leptonic sector of the electroweak theory. We 
shall use the effective four-Fermi low-energy description of that theory and we shall perform 
all our computations within the one loop (mean-field) approximation. Corrections to this 
picture from the inclusion of effects not captured by the low-energy effective theory, or 
from fluctuations about the mean-field approximation, must await further investigations. 
We shall discuss possible consequences of extensions of the dynamics beyond the standard 
model below. 

To set the stage, let us consider an action with a generic four Fermi coupling: 

s = ^VPlA^p - ij a Al^l] + MafaS^llp&tys ■ 

Here the index a summarizes all the attributes of ip: space-time dependence, gauge trans- 
formation properties, flavor, spin and whatever else. 

Because we are interested in pairings of the form (ipip) and (^'^'), we assume that M. 
has a Fierz-Bogoliubov (FB) decomposition: 

M a frs = 52v\Q { $Qps ) - 

A 

Here rjx = ±1, and Q^f = —Q$ because of the Fermi statistics of ip. Then we define 
auxiliary fields B^ x \ and add to C the term 

- E^(S (A)t - QftW(B (A) + Qfs%4>s) ■ 

A 

The original Lagrangian is recovered upon path-integrating on and . With the 
addition of this extra piece, the terms in £ that are quartic in tp and ^ cancel, leaving us 
with 
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S = -l^A*p - ^A T ^\ -J2V\B KX) B^ 

2 A 

+ V W + i)BH , 

A 

4* = -E^ t(A) g; ( 5 A) . 

A 

We obtain an effective action T e ff[B^ x \ B^] by integrating over ip and ip^. The result is: 
Te// = - E i7AS (A)t S (A) - + AA-'B^)-^} 

A 2 

where we have dropped some terms independent of B and B\ The mean field approximation 
that we shall adopt consists in demanding that T be stationary with respect to variations 
of B and £+: 

^ = ^ = 

These are the gap equations whose non-trivial solutions, if any determine whether a 
condensate with the quantum numbers of B^ and B^ can form. Because the effective 
theory we are dealing with is non-renormalizable, it is necessary to regularize these gap 
equations, and the size of the gap depends, apparently, on the value of the cutoff. Thus 
the gap equations in the form we derive them are not a good guide to the size of the gap. 
However, they allow us to determine which channels are attractive and hence in which 
channels one might expect a condensate to form. We shall see explicitly how this works in 
particular cases below. 

The simplest case is only one flavor of neutrino interacting with itself via the neutral 
current. Since vector exchange produces repulsion among like particles, we intuitively expect 
the v — v channel to be repulsive in this case. This is borne out by our explicit computation, 
which we now describe. 

The starting point is the four-Fermi interaction 

c int = -G^^yv 

where 75^ = — ip. The sign is dictated by the standard model. If we use 2-component 
fermions, and make use of the identity J2 a <T ai3 (T °5 = 2[ea 7 e<5/3 + ^apd^s} we see that 
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C int = 2G 2 (^ Jf a e a ^l f )(^/3e / 35^s) • (1) 

Making contact with our previous notation, we see that there is only one term in the 
FB decomposition with Q ai = y/2Ge ai , and rj — — 1. We remark that the fact that i] = — 1 
already tells us that the interaction is repulsive at the tree level; it is then to be expected 
that the gap equation will not have a non-trivial solution. 

For the kinetic part of the action, again using 2-component notation, we have A = 
ijfi — iff ■ V — A*) where ji is the chemical potential. Hence 



A -i = f d 4 P r ipo-v) -p-ff 

J (2vr) 4 (p - fi + iesgnpo) 2 - p 2 

and 



£ -ip-(x-y) 



[(P0+/i) ~p-ff T ] c -ip.( X -y) 

(27r) 4 (po + /j, + iesgnpo) 2 - p 2 



Note that the ie prescription has been introduced in the appropriate manner to take 
account of the role of \x as the chemical potential. One then proceeds to construct 
X = 4A~ 1 B(A T )~ 1 B' ( which, under the assumption (appropriate for a vacuum solution) 
that B and B^ are constants, can be written 

d P _-ip-(x-y) 



(2tt)< 

where 



+8G 2 B^B 

Hp) 



Po — (p-ff — AO 2 + it 

The gap equation will involve Tr(l + X)~ 1 X. Doing the po integral in this expression 
provides a factor of i that cancels the explicit i appearing in T e ff. The remaining integral 
over p is ultraviolet divergent. We regularize this by imposing a cutoff A on the magnitude 
of p. This leads to the unrenormalized gap equation, written in terms of M = B^ B 

-G 2 f A 1 
dpp 2 



7T 2 •'-a "* yj(p- fi) 2 + 8MG 2 

Clearly there is no solution (whereas there would have been a solution if we had had rj — 1; 
this would have changed the sign of the 1-h s). Note that for M — 0, there is a logarithmic 
divergence coming from p — /i. This is the singularity at the Fermi surface mentioned above. 
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Including an arbitrary number N of flavors turns out to be fairly straightforward assum- 
ing no flavor mixing. One begins with the analog of eqn. (1): 



-■int 



2GM (J Wf ) )(<^' ) ) (2) 



where i, j are flavor indices summed from 1 to N, and performs a separate FB transformation 
on the flavor indices using the identity 

2 Af2 ~ 1 

iV a=l 

where the A's are the generators of the fundamental representation of SU(N), normalized 
such that Tr X a X b = 25 ab . 

Following an analysis paralleling the one-flavor case£] one finds 

xU* -v) = -Ma J ^ T ^y~ lv ' {x ~ v) 



where 



4G 

T 



Po — {p-a — jj) 2 + ie 



and 



Mij = £ B™ A^Ag> . 

A,B 

Here the summation runs over the symmetric A's, including the unit matrix suitably nor- 
malized. It is important to note that M is a positive matrix: M = KK\ K = J2a B {A) X {A) . 
This expression for X leads to the following set of gap equations: 

B W = ^-B^ [ A dpp 2 tr[ , 1 . 



Multiply by B (A)t and sum on A. This produces 



■fin this analysis, we assume a common chemical potential for all the flavors. What may happen 
if this is not the case will be discussed briefly below. 
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V B ^ = / A dp p 2 tr[ . 1 M] 

a 27r • / - A ^/(p - /i) 2 + 4G 2 M 

whose only solution is = for all A. 

Thus, if we confine ourselves to neutrinos alone, and to the dynamics of the standard 
model, we find no possibility of neutrino pairing. Among the ways to avoid this conclusion 
are (a) extend the dynamics beyond the standard model (we shall discuss this possibility 
in the conclusions); (b) enlarge the dynamics to include the charged leptons (and possibly 
also the quarks). We have performed an analysis in which we consider not only neutrinos 
themselves but also electrons circulating in the loop. This generates an additional term in 
T e ff, and hence an additional contribution to the gap equation for the neutrino condensate. 
It does not, however, alter the result that there is no solution to the gap equation; (c) 
Finally, we can consider condensates that are composed not of neutrinos alone, but that pair 
neutrinos with charged leptons. Of course, since these condensates would be charged, their 
phenomenological consequences would be much more drastic than those of purely neutrino 
condensates. 

In any event, we begin with the Lagrange density C = £ + C int , where 
C = e(i ij- m)e + v e % ij v e + v^i ij v^- ft e e*e 

-g 2 - l. 

tint = ^2-K7 M (l - 7 5 )ee7^(l - 7 5 )z/ e - -^ e 7 M (l - ^)v e e^(l - 7 5 )e 

+ 2sm 2 6»H/e7 M ez/ e 7 M (l - 7 5 )z/ e - ^iyy^l - 75)^7^(1 - 7s)e 
+ 2sin 2 e w e^eiy ll -f ll (l - 75)^] . 

We have kept only the electron, as by far the lightest charged lepton. We have exhibited 
its couplings to u e and z/^; there should also be z/ T , which we ignore for simplicity. It enters 
in an identical manner to v^. 

Now we must implement the FB transformation. It is convenient first to make a Fierz 
transformation on the charged-current terms, to get them in a form resembling the neutral 
current terms, and then to make a FB transformation on the result. We find 

g 2 



w 



+ 4(1 + 2sm 2 6W)^ e 7°7 2 ez/ e 7°7 2 e 
+ 4(-l + 2sm 2 Mz7,7 7 2 e^7 7 2 e] 
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Note that 7°7 2 is the charge conjugation matrix, and hence expressions of the form z/7°7 2 e 
are proper Lorentz scalars. If we assume that the condensate will be a Lorentz scalar (this is 
only for convenience: there is no particular reason to expect that these condensates should 
respect Lorentz symmetry) then we limit ourselves to the last 2 terms. Note that they are 
of opposite sign. It turns out that the term involving v e is "bad" (77 = —1) whereas the 
one involving is "good" [r\ = +1). Hence we concentrate only on the last term, and 
introduce an auxiliary field B for it. We let k 2 = — 2sin 2 6 w ), and for notational 

ZM w 

convenience we introduce a doublet ip = (e, z/ M ) (a peculiar object from the point of view 
of the underlying standard model), so that, for example, mee = y^(l + r^)ip and p, e e^e + 
^v v \ v ii = i'K^i +^2^3)^. In our earlier notation, we have A = "f°(ip- y (1 + t 3 )) — Hi — fi 2 r 3 
and £ = _£b(i_ 75 ) 7 2 7 o Ti . 

Following the same analytic path as before, we arrive at a somewhat more complicated 
gap equation: 

l B*B 



2^tn IK B [ B r°° , 2 f°° , , 

kB , B = — / dp p / dp g(p ,p) 

7T 5 J-00 J-oo 



where 

G(po,p) 



p - fl+ + p 



[(Po - V+) 2 -P 2 - ™ 2 } [Po +H-+p]- 4k 4 BW[ Po -fi + +p] 



Here /i± — /ii ± /i2, and the integration over poles on the real po axis is by the prescription 
Po -> Po + iesgnpo. 

It is not possible to make much analytic headway with the expression in the general case. 
But if we set A 2 = 4k 4 B^B = 0, the expression becomes tractable, and we can then address 
two questions: (i) does the integral have the correct sign to permit a solution of the gap 
equation? and (ii) is there an infrared singularity, which would be evidence of an instability 
that could be cured by setting A 2 7^ 0? Let 00 = \/p 2 + m 2 > 0. Then the p integral can 
be done, yielding 



K 



2 rA . Of,.,* 



9(u 2 - ii 2 + ) (u - p)9(-p — //-) (UJ + p)9(p + fX-) - 



n 2 J-a uj u — p — 2fii uj + p + 2/Ji 

(u + p + 2fii)(u - p- 2/ii) 

Because of the ^-functions, each term in the integrand is non-negative, thereby answering 
the first equation in the affirmative. 
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Further, one sees that singularities can occur only at the Fermi momentum, p 2 F = n 2 _, and 
then only if the condition ji 2 + —m 2 = ji 2 _ is met. Since the electron density is proportional to 
[fi 2 + — m 2 fl 2 while the neutrino density is proportional to p?_ [9], this amounts to equating 
the Fermi momenta of the members of the pair. Phenomenologically, therefore, it is unlikely 
that this type of condensate could occur, except maybe in the early universe when larger 
background densities of both electrons and neutrinos were present. 

We turn now to a consideration of the size of the condensate. A simple BCS-\ike estimate 
gives 

A ~ p F e~^ . 

If the G 2 in eq. (1) or (2) is of order Gf, then, given the allowed range of pf, the 
exponential suppression makes A very small. One is led, therefore, to suggest the existence 
of a new interaction, acting only on neutrinos, for which the effective G 2 would have a scale 
of approximately 1 eV instead of the 200 GeV characteristic of Gp. An interaction of the 
same form as eqs. (1) or (2) would serve nicely, provided only we change the sign, thereby 
generating an attractive channel. The condensate would then produce Majorana neutrino 
masses of the form m v ~ A = G 2 (uu), and possibly a contribution to the cosmological 
constant A ~ G 2 | (vis) \ 2 . Since both G and the fermi momentum p F are of order 1 eV, one 
obtains neutrino masses and a cosmological constant that are likewise of this order. 

Furthermore, if we generalize our earlier analysis and allow the chemical potentials for 
the different neutrino species to vary, the condensates could depend non-trivially on flavor, 
perhaps leading to an interesting spectrum of neutrino masses and mixings. 

Our conclusions can be enumerated as follows: 

(i) There is no attractive channel in the purely neutrino sector of the standard model; 

(ii) The addition of charged leptons leads to attraction in the flavor off-diagonal channels, 
but a pairing instability occurs only if the Fermi momenta of the neutrino and the charged 
leptons are equal; 

(hi) In the case of neutrino-charged lepton pairing, there may also be the possibility of 
condensation in a Lorentz non-invariant channel. We have not looked at this in detail; 

(iv) If a new interaction exists among neutrinos with characteristic scale 1 eV, neutrino 
condensates could form with the right size to generate an interesting spectrum of masses 
and mixings, as well as an appropriate contribution to the cosmological constant. This 
possibility is currently under active investigation. 
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